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ABSTRACT 

Aims. We study the recently discovered twisting motion of bright penumbral filaments with the aim of constraining their geometry 
and the associated magnetic field. 

Methods. A large sunspot located 40° from disk center was observed at high resolution with the 1-m Swedish Solar Telescope. 
Inversions of multi-wavelength polarimetric data and speckle reconstructed time series of continuum images were used to determine 
proper motions, as well as the velocity and magnetic structure in penumbral filaments. 

Results. The continuum movie reveals apparent lateral motions of bright and dark structures inside bright filaments oriented parallel 
to the limb, confirming recent Hinode results. In these filaments we measure upflows of » 1.1 km/s on their limbward side and weak 
downflows on their centerward side. The magnetic field in them is significantly weaker and more horizontal than in the adjacent dark 
filaments. 

Conclusions. The data indicate the presence of vigorous convective rolls in filaments with a nearly horizontal magnetic field. These 
are separated by filaments harbouring stronger, more vertical fields. Because of reduced gas pressure, we see deeper into the latter. 
When observed near the limb, the disk-centerward side of the horizontal-field filaments appear bright due to the hot wall effect known 
from faculae. We estimate that the convective rolls transport most of the energy needed to explain the penumbral radiative flux. 

Key words. Sunspots - Sun: magnetic fields - Sun: photosphere — Techniques: spectroscopic - Techniques: polarimetric - 
Techniques: high angular resolution 



1. Introduction 

The discovery of twisting motions of penumbral filaments seen 
in time series of filtergrams by Ichimoto et al. (2007) ranks 
among the most striking recent discoveries in solar physics. In 
sunspots observed at a heliocentric angle 6 of 40 - 50° filaments 
lying roughly parallel to the limb display a twisting motion di- 
rected towards solar disk center. The nature and origin of this ap- 
parent motion is still unclear and requires further study. An im- 
portant step in this direction is to determine how these twisting 
filaments fit into the complex magnetic structure of the penum- 
bra. This is the main aim of the present letter, along with the con- 
firmation of the Hinode-based results of Ichimoto et al. (2007) 
with the 1-m Swedish Solar Telescope, which allows higher spa- 
tial resolution images to be obtained. 

2. Observations and data reduction 

Using the 1-m Swedish Solar Telescope (SST), we observed a 
mature sunspot in active region NOAA 10904 on August 13, 
2006 during a period of good to excellent seeing. The cen- 
ter of the field of view (FOV) was located at fi = 0.76 (he- 
liocentric angle 9 = 40.5°). By means of a dichroic beam- 
splitter, the sunlight was split into red (A > 500nm) and blue 
(A < 500nm) beams. In the blue channel, G-band (Ao = 430.5nm 
interference filter, FWHM=1.3 nm), and blue continuum (Ao = 
436.4nm, FWHM=1.1 nm), images were simultaneously col- 
lected on Kodak Megaplus 1.6 CCD cameras (1536x1024 pix- 
els, of 9 fim pixel-size). The exposure time was set to 13 ms. 
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These images were binned into packages of 60 frames and pro- 
cessed using speckle image reconstruction techniques (Weigelt 
1977; Pehlemann & von der Liihe 1989; de Boer 1996). By 
applying this procedure, we obtained a two-hour long time se- 
quence of near-diffraction-limited filtergrams at a cadence of 
19 s between individual frames. 

In the red channel three temporally synchronized Sarnoff 
CAM 1 Ml 00 cameras were employed, with exposure times 
of 4.5 ms. Two of them collected broad-band continuum 
frames near A = 630.25 nm. The third camera registered 
spectro-polarimetric data consisting of scans through the Fe I 
(A = 630.25 nm) line at 6 wavelength positions (A - Ao = 
[-150,-75,0,75,150,250] mA) made with the Solar Optical 
Universal Polarimeter (SOUP) filter. The last position in general 
samples the continuum. Two Liquid Crystal Variable Retarders 
(LCVRs) were used to modulate the light beam. The complete 
scan across the iron line lasted about 123 s. After calibration 
procedures, the full Stokes vector at each point of the FOV was 
retrieved. For this we used a telescope model and demodulation 
matrices of the optical setup measured with dedicated calibration 
optics (Selbing 2005). 

We determined the magnetic field vector and the line-of- 
sight (LOS) velocity by inverting the measured Stokes profiles 
(composed of only 6 wavelength points) using the HeLIx inver- 
sion code (Lagg et al. 2004). An atmosphere composed of two 
Milne-Eddington atmospheric components, one magnetic, the 
other one representing unpolarized stray light, was employed. 
Extensive tests based on inverting synthetic spectra computed in 
3-D radiation MHD models have shown that 6 wavelength points 
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Fig. 1. Temporal evolution of an individual filament in blue con- 
tinuum images. The size of the FOV is 1.5x3.5 arcsec and the 
cadence is Af = 76 s. White arrows mark the position of a bright 
structure that moves from left to right. The solar limb is located 
to the left, the umbra to the bottom. 

are sufficient for a reliable inversion. Q The inversions provide 
the magnetic field vector: (field strength |B, m |, inclination y, m ., 
azimuth 0; nv ) and the line-of-sight (LOS) velocity, vlos, (which 
was calibrated to be m/s in the darkest part of the umbra). We 
applied an algorithm based on minimization of the gradients of 
the azimuth over the FOV to correct for the 180° ambiguity in 
</>,„,,. Additionally, B, m , was constrained to point away from the 
sunspot umbra due to its positive polarity (confirmed by MDI 
magnetograms). If even after this algorithm for a given pixel i, 
> 90° (A0 /m , = \(f> invJ - 0, m ,,+i|), with respect to the (f> inv 
of 8 neighbouring pixels, then we considered the results of the 
inversion of pixel i to be incorrect. The values in this pixel were 
then replaced by the median of the surrounding pixels. To anal- 
yse the geometry of the magnetic field better, a rotation matrix 
was applied to B, m , to get the magnetic field vector B - (B, y, </>), 
in local solar coordinates. 

3. Results 

Let us first consider the blue continuum time sequence, which 
has the best spatial resolution of our data. The most remarkable 
of the rapid intensity variations within penumbral filaments are 
continuous lateral motions of dark and bright stripes across sev- 
eral bright filaments, giving a visual impression of a twisting 
motion or a rotation around their axes. To illustrate the temporal 
evolution of one such filament located in the central penumbra, 
we plot in Fig. 1 a 228 s sequence of 4 blue continuum images 
(1.5x3.5 arcsec in size), with a cadence of 76 s. The filament 
makes an angle of 5 — -21° to the nearest limb (a negative num- 
ber implies that the outer part points towards disk center). We 
filtered out the typical jitter of the images, which comes from 
image distortions due to residual seeing effects, using a subsonic 
filter that cuts out all motions faster than 4 km/s. In Fig. 1 an 
elongated bright stripe moves from the left side to the right of the 
filament, i.e. towards disk center. Another dark and bright struc- 
tures in this and other filaments move in the same direction. The 
average horizontal velocity of such motions is « 1 km/s. These 
lateral motions are always directed towards the disk-center side 
of the filaments, regardless of their position in the penumbra 
(with respect to the umbra). These results agree with those ob- 
tained by Ichimoto et al. (2007). Such motions can be identified 

1 These tests were carried out by A. Lagg et al. (unpublished) in the 
course of studying the design requirements of the Visible-light Imager 
& Magnetograph for ESA's Solar Orbiter Mission. 



over a large part of the penumbra. They are most prominent in 
filaments oriented nearly parallel to the nearest limb (6 = 0° or 
180°). As the azimuth of the filaments to this direction increases, 
the lateral motions become less visible. The effect vanishes at 
angles -140° < 6 < -40° and 40° < 6 < 140°. 

Next, we consider the results of inversions of the Stokes pa- 
rameters. In Fig. 2 we plot maps of normalized continuum in- 
tensity near A = 630.25nm, I c lh (where Iq means the contin- 
uum intensity averaged over the quiet Sun present in the same 
frame), y, <p, and Vios> of the same region as in Fig. 1. The 
magnetic vector is given in solar coordinates. The dashed cen- 
tral contour outlines the bright part of the filament {I c > 1 .0) 
as is visible at the chosen viewing geometry. The maps of |£|, 
y, </>, and \'ios display an offset relative to the intensity map. 
Partly this is because these are deduced from the A = 630.25nm 
line whose flanks are formed typically 150 km above t c — 1 
in the penumbra, as we estimated from contribution functions 
convolved with the filter profiles. To compensate for this par- 
allax, we need to shift the continuum contour by a; 100 km to 
the left in the 4 right panels in Fig. 2, which results in the solid 
contours. Even after such a shift, a significant offset remains, 
which has a physical origin. Thus, the limb side of the filament 
shows blue-shifts, whereas its center side displays weaker red- 
shifts. The center-ward side of the filament has enhanced mag- 
netic field strength and greatly reduced zenith angle (i.e. a more 
vertical field), whereas an the limbward side the magnetic field 
is weaker and is almost horizontal. The azimuth at the limb side 
of the filaments is higher than at the disk-center side. Note that 
= corresponds to a direction parallel to the jt-axis and points 
towards the right, so that at the disk-center side the projection of 
the field is roughly parallel to the filament, while on the limb- 
ward side it points « 20° in the direction of the limb. 

Figure 3 displays a horizontal cut, made along the dashed 
line in Fig. 1 (frame 4) and in Fig. 2, of physical parameters. 
Figure 3 a illustrates the evolution of the blue continuum inten- 
sity along the cut. Lower curves were recorded later at a cadence 
of 38 s. Fig. 3b shows the variations of vlos- The velocity pro- 
file shows an upflow (around —1.13 km/s) on the limb side and 
a weak downflow (around +0.09 km/s) on the center side of the 
filament. This asymmetry between up- and downflows is also 
displayed by other studied filaments. The velocity profiles pre- 
sented by Ichimoto et al. (2007) are also strongly asymmetric, 
i.e. display a strong blue shift and a weak red shift. This suggests 
that the asymmetry is an intrinsic property of the filaments, al- 
though we cannot rule out that it is an artifact due to crosstalk 
from the Evershed flow, a small component which may point in 
the direction of the LOS. 

The continuum intensity (near/l = 630.25nm) too has been 
horizontally shifted by 100 km to allow a parallax-free compari- 
son with velocity and magnetic field vector. The different spatial 
resolution and sensitivity to contrast in the blue and red chan- 
nels, as well as the fact that they were recorded at different times, 
may be responsible for the rather distinct continuum profiles in 
Figs. 3a and 3b. The center of the bright filament nearly coin- 
cides with the highest blueshift. These results are consistent with 
the findings of Ichimoto et al. (2007). Figures 3c,d,e display the 
components of the magnetic vector |H|, y, and </>, respectively. 
The magnetic field strength reaches its maximum about 0.3 arc- 
sec towards disk center from the axis of the bright filament. At 
this point B is also most vertical (y =28°) and its azimuth is 
around 86°, which corresponds to an orientation of the surface 
projection of the magnetic vector, B SU rf. roughly parallel to the 
filament. On the limb side of the bright filament and in the dark 
region immediately bordering it, \B\ is nearly a factor of 2 lower 
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Fig. 2. From left to right: continuum intensity near 630.25 nm normalized to its value averaged over the nearest quiet Sun area, 
magnetic field strength \B\ [G], zenith angle, y [°], azimuth angle in the solar surface plane, <p [°], and the line-of-sight (LOS) velocity, 
vlos [m/s], respectively. Dashed contours localize the bright penumbral structures (I c > 1.0). Solid contours have been shifted to 
roughly compensate for the different heights to which continuum and the other physical parameters refer. (f> = corresponds to a 
direction parallel to the X-axis and points towards the right. 



and B becomes almost horizontal (70° < y < 90°). Here the az- 
imuth is skewed by 15° < (f> < 25° towards the limb relative to 
the axis of the filament. 



4. Discussion 

In Fig. 4 we present a sketch of the geometry of the penum- 
bral fine structure, which is consistent with the results of the 
present letter, as well as those of Ichimoto et al. (2007), Borrero 
et al. (2008), and Borrero & Solanki (2008) based on Hinode 
data. The r = 1 level is relevant for the line wings. In the 
sketched configuration stronger, more vertical fields (spines) are 
associated with less dense gas, leading to a lowered T = 1 
level there. Interspersed between them are filaments of nearly 
horizontal, weaker magnetic field pointing into the page (inter- 
spines). Important is that the filaments carrying the Evershed 
flow have a sizable horizontal field (Borrero & Solanki 2008), 
since the flowing gas is observed to be magnetized (Solanki et al. 
1994; Bellot Rubio et al. 2004). It is within these filaments that 
the transverse motions are seen, which are similar to overturn- 
ing convection. A nearly horizontal magnetic field means that 
convection takes place in the form of convective rolls, as first 
proposed by Danielson (1961). Due to the raised r = 1 level in 
these filaments (consistent with the findings of Rimmele 2008) 
the bright part of the filament, which is contoured in Fig. 2, is the 
penumbral counterpart of the hot wall found in faculae and pores 
(Spruit 1976; Keller et al. 2004; Carlsson et al. 2004; Cameron 
et al. 2007). This geometry explains why only lateral motions 
directed towards disk center are seen and why the observed up- 
flow velocity is higher than the downflow velocity, as can be 
seen by considering the directions of the flow arrows relative to 
the line-of-sight in Fig. 4. From these hot walls radiation flows 
into the neighbouring spines. The parallax effect, discussed in 
Sect. 3, makes it obvious that the most inclined magnetic field, 
which is located above the center of the filament, was detected 
at the limbward border of the bright filament. The contrast be- 
tween spine and interspine is expected to depend on the rela- 
tive temperature contrast, width of the spine and the difference 
in evacuation between spines and interspines. The values plot- 
ted in Figs. 2 and 3 are weighted averages along the LOS (the 
relatively few measured wavelength points do not allow a more 



detailed analysis). Thus, because a number of the slanted rays 
pass through both the inclined and the horizontal fields, the dif- 
ference between them appears less clear in Figs. 2 and 3 than 
it may be in reality. Whether our measurements contradict the 
findings of Jurcak & Bellot Rubio (2008), who report that bright 
penumbral filaments show the more vertical fields and weaker 
flows, cannot be judged, in part because of the different view- 
ing geometry. They have investigated the limbward part of the 
outer penumbra, whereas our investigation addresses filaments 
oriented parallel to the limb, in the mid penumbra. 

The energy transport by the convective rolls can be ex- 
pressed as F CI = p ■ u ■ v, (neglecting the enthalpy) where 
p = 2 x 10~ 7 g crrT 3 is the mass density, u = 3 x 10 12 erg g -1 
is the heat deposited by 1 g of gas as it cools from 12000 K 
to 5000 K (Schlichenmaier et al. 1999), and v = 1 km s _I is 
the upflow velocity observed in the convective rolls. Entering 
these numbers, we obtain F a w 6 x 10 10 erg cirT 2 s . This is 
larger than the radiative flux emitted by the penumbra, F pen ss 
4.7 x 10 10 ergcirT 2 s _1 = 0.75 ■ Fq. However, we must keep 
in mind that upflows fill at the most half of the surface area of 
the penumbra; e.g., the spines must be heated radiatively. It is 
reasonable to assume them to emit the same flux as the umbra 
in the absence of any lateral heat influx from the interspines, 
F U mb = 0.2 • Fq. For a penumbra half covered by upflows and 
half by gas at umbral temperatures, we obtain the following re- 
lationship that must be fulfilled: j ■ {p ■ u • v + F um b) = F pen . This 
gives a requirement on p ■ u ■ v of « 8 x 10 10 ergcirT 2 s , 
which is only a factor of 1 .3 larger than the estimate obtained 
from the convective rolls. An underestimation of the partly un- 
resolved velocities of the roll convection can easily account for 
this factor. The observed roll convection transports far more en- 
ergy flux than interchange convection or the Evershed flow along 
flux tubes, as estimated by Schlichenmaier & Solanki (2003). We 
therefore propose that the observed convective rolls are the main 
form of energy transport in the immediate subsurface layers of 
the penumbra, carrying most of the energy required to maintain 
the radiative output of the penumbra. 

The configuration of field and flows presented here combines 
aspects of both (i) the uncombed penumbra model of Solanki & 
Montavon ( 1 993) and (ii) Schlichenmaier et al. ( 1 998a,b) and the 
gappy penumbra model of Spruit & Scharmer (2006), Scharmer 
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& Spruit (2006), cf. Scharmer et al. (2008). The geometry of 
magnetoconvection in the umbral dots models by Schiissler & 
Vogler (2006), cf. Riethmuller et al. (2008), displays qualitative 
similarities to our model. Observations by Rimmele (2008) of 
a sunspot close to disk center reveal that the velocity pattern of 
dark-cored filaments at lower layers is similar to what he sees 
in the light bridges, i.e. the dark lanes are correlated with up- 
flows and with downflows to the sides. This also agrees with our 
results. 
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Fig. 3. From top to bottom: a) horizontal cuts of blue contin- 
uum intensity, b) LOS velocity (solid) and red continuum inten- 
sity (dashed), c) strength, d) zenith angle, and e) azimuth of the 
magnetic field inside the filament. The solar limb is located to 
the left. The continua profiles have been shiftet to the left by 
0.1 arcsec. In panel a) intensity curves recorded later are off- 
set downwards (cadence = 38s). Diamonds mark the location 
of a bright structure drifting across the filament. Panels b-e dis- 
play the results of data recorded between 09:16:50 - 09:18:30. 
The vertical dashed-triple-dotted line corresponds to the loca- 
tion, averaged over the 09: 16:50 - 09: 18:30 time interval, of red 
continuum intensity maximum of the bright filament. 
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Fig. 4. Schematic illustration of the field and flow structure in 
penumbral filaments. Plotted is a vertical cut placed perpendicu- 
larly to the filaments (white semi circular shapes at the bottom of 
figure). The thick dashed-dotted line designates the r = 1 level 
for the plotted view direction. Vertical arrows denote magnetic 
field lines and arrows inside the lightened filaments denote con- 
vective flows. Circled crosses represent magnetic fields inside 
the filaments oriented along their axes (which also correspond to 
the direction of the Evershed flow). 



